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Abstract 

We show that the free energy of the three-state T%{t q ) model can be ex- 
pressed as products of Jacobi elliptic functions, the arguments being those of 
an hyperelliptic parametrization of the associated chiral Potts model. This is 
the first application of such a parametrization to the A-state chiral Potts free 
energy problem for N > 2. 
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1. Introduction 

In the field of solvable models in statistical mechanics, the A-state two- 
dimensional chiral Potts model has proved particularly challenging. This is 
despite the fact that for A = 2 it reduces to the Ising model, the free energy 
of which was calculated by Onsager (1944). 

The model was developed by Howes, Kadanoff and den Nijs (1983), von 
Gehlen and Rittenberg (1985), Au-Yang et al (1987) and McCoy, Perk and 
Tang (1987). It was first fully defined as a general A-state lattice model 
by Baxter, Perk & Au-Yang (1988), who showed that it satisfied the star- 
triangle relations. These relations define "rapidity" variables p, q such that 
the Boltzmann weights are functions W pq of p and q. 

Many of the citations herein are to earlier papers by the author: in these 
we shall abbreviate "Baxter" to "B". 

With previously solved models, such as the hard-hexagon model (B 1980), 
the calculation of the free energy was straightforward once one had obtained 
the star-triangle relation and rapidities for that model. Typically, it turned 
out that there was a transformation in terms of Jacobi elliptic functions such 
that Wpq depended on p, q only via their difference u = q — p. The free energy 
is — ksT log n pq , where K pq is the partition function per site. This K pq must 
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also be a function k(u) of u, and there were always "rotation" and "inversion" 
relations (B 1982a, 1982b), of the form 

k(u) = k(\ — u) 

k(u) k{—u) = r{u) , (1-1) 

where A is a positive real constant (the "crossing parameter") and r(u) a 
known meromorphic function. For u real and < u < A the Boltzmann 
weights are real and positive, and one can develop low-temperature series ex- 
pansions in the usual way that indicate that k(u) is analytic, non-zero and 
bounded in the vertical strip < 5ft(u) < A. The equations (|l.lj) then deter- 
mine k(u) and one can solve them by Fourier transforms or other methods. 

The N > 2 chiral Potts model, however, does not have the "rapidity dif- 
ference property" , i.e. there is no transformation that takes W pq to a function 
only oiq—p. This makes the calculation of n pq much more difficult, and it was 
not till 1990 that an explicit result was obtained (B 1990, 1991a). The calcu- 
lation of the spontaneous magnetizations (order parameters) is even harder, 
and was not accomplished until 2005 (B 2005a, 2005b). 

The calculation of n pq proceeds in two stages. First one calculates the 
partition function per site T2(p,q) of an associated u T2(t q ) model", which is 
closely related to the superintegrable case of the chiral Potts model. Then 
one uses this result to calculate tz pq . 

In the calculation of the order parameters, we used the fact that certain 
functions G p y p (r), S(p) are similar to T2(p,q). In fact they are ratios of 
special cases of T2(p, q). 

Although the main task, namely the calculation of the free energy and or- 
der parameters of the chiral Potts model, has been accomplished, it remains 
disappointing that one has no elegant parametrization in terms of elliptic func- 
tions, as one has for models with the difference property. These parametriza- 
tions explicitly exhibit the poles and zeros of n pq on an extended p, q Riemann 
surface. There is a parametrization of the rapidities and Boltzmann weights of 
the A^-state chiral Potts model in terms of hyperelliptic functions (B 1991b), 
but these have N — 1 arguments that are related to one another in a com- 
plicated way, and until recently they have not been found to be particularly 
useful. 

However, we have now found that for = 3 the function S(p) mentioned 
above can be expressed as a ratio of generalized elliptic functions of these 
arguments (B 2006). The obvious question is whether T^ip-, q) can be similarly 
expressed. We show here that the answer is yes. The functions that occur are 
the same as that appear in other solvable models, notably the Ising model. 

2. The function T2(p, q) 

For the chiral Potts model, the rapidity p can be thought of as the set of 
variables p = {x p ,y p , fi p ,t p }, related to one another by 

t p = x p y p , x p + y p = k{\ + x p y p ) , 
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(2.1) 

= 1 k I , ^Up == 1 ^ l^p 

Here k, k' are real positive constants, satisfying 

k 2 + k' 2 = 1 . (2.2) 

In terms of the a p ,bp,c p ,dp of (Baxter, Perk & Au-Yang 1988; B 1991b), 
Xp — dp/dp, Up — bp/ Cp, Hp — dp/ Cp. 

There are various automorphisms or maps that take one set {x p , y p , n p ,t p } 
to another set satisfying the same relations ()2.1j) . Four that we shall use are: 

R ■ {xR p ,yRp,fJ,R P ,t R p} = {y p ,u;xp,l/fip,ujtp} , 

U : {xup,yup,fJ-Up,tu P } = {x^ 1 ,y~ l ,u~ l/2 Xpiip/yp,t~ 1 } , 

V : {xvp,yv P ,^Vp,tVp} = {x p ,u;yp,Hp,ujtp} , (2.3) 

M : {xm P ,VM P , ^Mp,tMp} = {Xp,y p ,LOHp,tp} , 



where 
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They satisfy 



RV~ l R = V , MRU = R 



U 2 M = V N = M N = 1 . (2.4) 

Here U = RSV, S being the operator S used in (Baxter, Perk & Au-Yang 
1988, B 2006). 

We take [i v to be outside the unit circle, so 

IMpI > 1 • (2.5) 

Then we can specify x p uniquely by requiring that 

- vr/(2iV) < arg(xp) < ir/(2N) . (2.6) 

We regard Xp,y p ,/j,p as functions of t p . Then t p lies in a complex plane 
containing N branch cuts £>o, Si, ... , Bjv-i on the lines arg(i p ) = 0, 2tt/N, . . . , 
2ir(N — 1)/A r , as indicated in Fig. ^ while x p lies in a near-circular region 
round the point x p = 1 , as indicated schematically by the region TZq inside the 
dotted curve of Fig. ^ The variable y p can lie anywhere in the complex plane 
except in TZq and in N — 1 corresponding near-circular regions TZi, . . . , 1Zn~i 
round the other branch cuts. With these choices, we say that p lies in the 
"domain" V. 

It can be helpful to consider the low-temperature limit, when k' — > 0. The 
branch cuts Bq, B\, ... , Bn-i and the regions 1Zq,1Zi, . . . , Hn-i then shrink 
to the points 1, u>, . . . , lo n ~ x . If t p is held fixed, not at 1, u>, . . . , w^ -1 , then 
z P -> *p and /x^ = 0(l/k'). 

Define q similarly to p, also lying in P. Then from eq. 39 of (B 1991a) 
and eq. 54 of (B 2003a), the partition function per site T 2 {p, q) = T2(fi p , t g ) of 
the T2(t q ) model is given by 

lo g r 2 (p,q) = _jf ^_ el ^ j log [A(0) - wt,] d# , (2.7) 
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Figure 1: The cut t p -plane for N = 3, showing the three branch cuts B ,Bi, B 2 and 
the approximately circular region TZq in which x p lies when p G T>. 



where 



A(0) = [(l-2k'cos6 + k' )/k 2 ] 



2]1/N 



(2.8) 



We have omitted a factor y 2 from t 2 (p, q) , and to ensure that this formula is 
correct for p,q eP, we invert fj, q in (B 1991a, 2003a). 
This function T 2 (p,q) satisfies the relations 



and 



where 



T2(Vp,q) = T 2 (p,q) , 
N-l N-l 

II r 2 (p,V j q) = J] M» p ,iJt q ) 
i=o j=o 



a(p, g) 



a(p, g) 



AT _ „,N\( X N 



(2.9) 
(2.10) 

(2.11) 



3. The Riemann sheets ("domains") 
formed by analytic continuation 

We shall want to consider the analytic continuation of certain functions 
of t p onto other Riemann sheets, i.e. beyond the domain T>. We restrict 
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attention to functions that are meromorphic and single- valued in the cut plane 
of Figure n an d similarly for their analytic continuations. Obvious examples 
are x p ,y p and T2(t p ). They are therefore meromorphic and single- valued on 
their Riemann surfaces, but we need to know what these surfaces are. 

We start by considering the most general such surface. As a first step, allow 
\i v to move from outside the unit circle to inside. Then t p will cross one of 
the N branch cuts Bi in Figure Q moving onto another Riemann sheet, going 
back to its original value but now with y p in TZi. Since y p is thereby confined 
to the region near and surrounding oj 1 , we say that y p ~ oj 1 . Conversely, by 
y p ~ uj' l we mean that y p € TZi 

We say that p has moved into the domain T>i adjacent to V. There are N 
such domains Vq, V\, . . . , Djv_i. 

Now allow fip to become larger than one, so t p again crosses one of the 
TV branch cuts. Again we require that t p returns to its original value. If it 
crosses Bi, then it moves back to the original domain T>. However, if it crosses 
another cut Bj then x p moves into TZj-i, and we say that p is now in domain 

Proceeding in this way, we build up a Cayley tree of domains. For instance, 
the domain T>ijk is a third neighbour of T>, linked via the first neighbour T>i 
and the second-neighbour T>%j, as indicated in Figure EJ Here x p ~ 1 in V, 
y p ~ lo 1 in T>i, x p ~ w J in T>ij and y p ~ u> k in T)^. We reject moves that 
take p back to the domain immediately before the last, so j ^ and k ^ i. 
We refer to the sequence {i,j, k, . . .} that define any domain as a route. We 
can think of it as a sequence of points, all with the same value of t p , on the 
successive Reimann sheets or domains. 

The domains T>, Vij, T)^^... with an even number of indices, have x ~ u , 
where £ is the last index. We refer to them as being of even parity and of type 
t. The domains T>i, T>^,... have y ~ u/ and are of odd parity and type I. 



V V>i T>ij Vijk 



Figure 2: A sequence of adjacent domains V, Vi, Vij,Vij k . 

The automorphism that takes a point p in T> to a point in X>j, respectively, 
is the mapping 

Ai = V 1 ' 1 RV~ l . (3.1) 

If p' = Ai p, then 

Xp> = to~ l y p , y p > = uj l x p , t p > = t p . (3.2) 
Because of (|2.4|) . Aj+jy = Ai, so there are N such automorphisms. 



We can use these maps to generate all the sheets in the full Cayley tree. 
Suppose we have a domain with route {i, j, k, . . .} and we apply the automor- 
phism A a to all points on the route. From IJM.2JI this will generate a new route 
{a, i — a, j + a, k — a, . . .}. For instance, if we apply the map A a to the route 
{m} from T> to T> m , we obtain the route {a,m — a} to the domain D ajrrir<x . 
Thus the map that takes T> to T>ij is AiAi + j. 

Iterating, we find that the map that takes T> to T>ijk... m n is 

A-iAi+jAj+k ■ ■ ■ A m + n . (3-3) 

We must have 

A 2 , = 1 , (3.4) 

since applying the same map twice merely returns p to the previous domain. 

Let us refer to the general Riemann surface we have just described as Q. 
It consists of infinitely many Riemann sheets, each sheet corresponding to a 
site on a Cayley tree, adjacent sheets corresponding to adjacent points on the 
tree. A Cayley tree is a huge graph: it contains no circuits and is infinitely 
dimensional, needing infinitely many integers to specify all its sites. 

Any given function will have a Riemann surface that can be obtained from 
Q by identifying certain sites with one another, thereby creating circuits and 
usually reducing the graph to one of finite dimensionality. 

From (JH2J), the maps Aq,A\,... , An leave t p unchanged. We shall often 
find it helpful to regard t p as a fixed complex number, the same in all domains, 
and to consider the corresponding values of x p , y p (and the hyperelliptic vari- 
ables z p , w p ) in the various domains. To within factors of u, the variables x p 
and y p will be the same as those for T> in even domains, while they will be 
interchanged on odd domains. 

(a) Analytic continuation ofr2(p,q) 

Consider T2(p,q) as a function of q (so replace p by q in the previous 
discussion of the domains). More specifically, think of it as a function of the 
complex variable t q . For q € T> (and p £ D), it is apparent from (|2.7|) that 
T 2(p,q) is an analytic function of t q except for the single branch cut B^-i-, 
being single- valued across the other N — 1 cuts. Let q' = Ai q, so q' G T>i, and 
define 

=1 if i = j, mod N, 
5(i,j) = else. (3.5) 

Then it follows from ((2~TUj) that 

TiM) = v m 5 ^ N ~^T 2 (p,q) (3.6) 

for i = 0, 1, . . . , N — 1, where 



a(p,q>) ( x N_ x N ){y N_ AT) 



l/v pq > (3.7) 



and = yq,y$ = x*. 
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Eqn. (|3.6|) defines the mapping Aj applied to the function T2(p,q) of q. 
Iterating, it follows that if q" = Ai Aj A^ q, then 

T 2 (p,q") = v M m T 2 (p,q) , (3.8) 

where m = 5(i, N - 1) - 5(j, N-l) + 5(k, N - 1). 

We can also keep q fixed in T> and consider the analytic continuation of 
T 2(p, q) as P moves from sheet to sheet. If p' = Aip, so p £ T>, p' € T>i, we can 
verify from (|2.7f) that 

T2(p',q) = (uj~H p - UJtq) 2 / r 2 (p, q) , t p i=tp . (3.9) 

Also, from (|2,lll) . 

a(p,q)a(p',q) = (t^-t^) 2 . (3.10) 

Eqn. (|3.9|) defines the Aj for the function T 2 {p, q) of p. Iterating, it follows 
that if p" = Ai Aj A/, p, then 

T 2 {p",q) = (u;-%-ut q ) 2 {u-H p -Lot q )- 2 (io- k t p -iot q ) 2 /T 2 {p,q) . (3.11) 

Note that for both T 2 (p,q) —> T 2 (p,q") and T 2 (p,q) —>■ T 2 (p",q), it is true 

4. Hyperelliptic parametrization for = 3 

Hereinafter we restrict our attention to the case iV = 3. A parametrization 
of {x p , y p , /ip, t p } was developed in previous papers, in terms of a "nome" x 
and two related parameters z p ,w p (B 1991b, 1993a, 1993b, 1998). The nome 
x is like k and k' in that it is a constant: it is not to be confused with the 
rapidity variable x p . 

However, in (B 2006) we showed that the function S p could not be ex- 
pressed as a single- valued function of these original parameters z p ,w p . This 
is because z p ,w p have the same values (for given t p ) in the domains T>o2i and 
D21I) whereas S p has different values therein. These domains are obtained by 
the maps AqA 2 Aq, A 2 AqA 2 , respectively, and indeed we see from (|3.8j) and 
(|3.9|) that these maps give different results for r 2 (p, q) for both the q and the 
p variables. Thus T 2 (p,q) cannot be a single-valued function, either of z q ,w q 
for fixed p, or of z p , w p for fixed q. 

The same problem occurs with the domains T>\ 10, D220 an d the corre- 
sponding maps A\A 2 A\, A 2 A\A 2 . 

The situation is not lost. We also showed in (B 2006) that there is an- 
other way of parametrizing k, x p ,y p , fJ> p ,t p . This alternative way preserves the 
property that the nome x is small at low temperatures (k' small). It can be 
obtained from the original parametrization by leaving unchanged and 

transforming k, k', x p , y p , fi p , t p according to the rule: 

k, k , Xp, yp, Hp, tp > k , ik jk, 1/ Xp, yp, lj ^ Xp^ip, yp/ Xp (4.1) 

(taking uj~ 1 / a = e~ m / 2N ). This mapping leaves the relations (|2.1|) . (|2.2() 
unchanged. 
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Doing this, eqn. (21) of (B 1993a) becomes 

oo / 1 _ 3n \ 12 

- k ' 2 = ^U(y^r) > ( 4 - 2 ) 



while the two equivalent relations (4.5), (4.6) of (B 1993b) remain unchanged: 
00 (1 - x 2n -V<u;)(l - x 2n - l w/z)(l - x 6n - 5 zw)(l - x 6 ™- 1 ^- 1 ™- 1 ) 



w 



- x (1 - x 2n - 2 z/u!)(l - x 2n u>/z)(l - x 6n ~ 2 zw)(l - x^-^z^w- 1 ) 



Z _yr (1 - X 2n " 2 »(l - X 2n W )(l ~ X 6n - 4 Z 2 / W ){l - X^W/Z 2 ) 

w ~ (1 - x 2n ~ 1 t(;)(l - 2; 2n - 1 /w)(l - x 6n - 5 w/z 2 )(l - x en ~ 1 z 2 /w) 

(4.3) 

These z,w are related to x p ,y p , fj, p ,t p by various elliptic-type equations 
that we shall give below, being the arguments (more precisely the exponentials 
of i times the arguments) of Jacobi elliptic functionsof nome x. Thus x, like 
k,k', is a constant, while z,w are two more rapidity variables, dependent on 
p. We shall them as z p , w p . 

First we introduce elliptic-type functions 

ur ^ ,/ur-U 2w v 2 ^ (l-uJX n - 1 z)(l-L0 2 x n /z) 
h(z) = 1 h(z 1 )=u 2 h(xz) = -u 2 [[ )- / ^j-L , 

(1 — UJ Z X n 1 Z)(l — LOX n jz) 

™ (l-x 3n ~ 2 /z)(l-x 3n ~ 1 z) 
~ l} l {l-x 3n - 2 z){l-x 3n ~ 1 /z) ' 

<l>{z) = l/^(z~ l ) = <p(x 3 z) = z^Mz) (4-4) 

so that h(l) = 0(1) = 1, and define two sets of parameters p,p in terms of 
z p , w p . 

Pi = z p , p 2 = -l/w p , P3 = -w p /z p , (4.5) 

Pj = Pj+l/pj-l , 

extending pj , pj to all j by 

Pj+3 = Pj , Pj+3 = Pj • (4.6) 

Thus 

p x = z p /w 2 p , p 2 = -w p /z p , p 3 = -z p w p . (4.7) 
We shall also need certain cube roots pj of x pj , choosing them so that 

Pj = ( x Pj) 1/3 = ^Pj+iPj-i ■ (4.8) 

Then after applying the mapping (|4.1|) . the equations (27), (32) of (B 1993a) 
become 

y p /x p = -u~ j h(pj) , (4.9) 
= ~u~ 3 (l>(xpj) = -uj-Jpjfoixpj) , (4.10) 

X p [l p 
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for j = 1,2,3. The second form of (|4.1U|) is to be preferred as it fixes the 
choice of the leading cube root factor in the function (f>(z), so fixing pj. 

Similarly, replacing p by q in the above equations, we define qj,~q~j,qj, and 
hence x q ,y q , fi q ,t q , in terms of two related variables z q ,w q . 

(a) The low-temperature limit 

At low temperatures x, k' are small. For p £ T> (and t p not close to a cube 
root of unity) we can choose z p , w p to tend to non-zero limits as x — > 0. Then 
(|4.3|) both give w p = z p + 1. Also, x p — > 1, so Q2.1J) and (|4.9|) determine x p , y p 
uniquely (with the same value for j = 1, 2, 3). Then /ip can be calculated from 
the last of the eqns. I|2.1[) . and jjL p itself from Q4.1UJI . Hence when x,k' are 
small 

k — 2 7x ; — 2^ "I - 1 ? — 1 ? 

2 

2/ P = £ , k'fi p = -3(1 -u; J y p )pj , 

w — 

for j = 1, 2, 3. Note that x p ,y p , z p ,w p , Pj,Pj are of order unity, but the pj are 
of order x 1 / 3 , while n p is of order a; -1 / 6 . 

(6) Mappings 

The effect of the automorphisms R, U, V on z p , w p is 

£Rp = -XWp , Z Up = -W p , Z Vp = -1/Wp , 

WRp = Wp/Zp , W Up = -Z p , W Vp = Zp/Wp . (4.11) 

The mappings U, V, M take a point p within T> to another point within T>: 

p' = Up: p'j = l/p 3 -j, Pj = P-j , Pj = u^p-j , 
P ' = Vp: p'j = p j+1 , p'j = p j+1 , p'j = p j+ i , 
p' = Mp: p'j = Pj , p'j = Pj , Pj = ufij ■ (4.12) 

and if p' = Rp, then p'j = u> 2 ~^ x l ~ 5 ^' 2 ^> /pj+±. 
Also, if p' = Aip, then 

Zp , = x^-^^/zp , W p , = X^/Wp 

p>. = x 6(i+j,2)S(i+j,l) / p . j (fi.y = u -i-3 x 1 -S(i+m/p j ; (4.13) 

for i = 0, 1, 2. We see that in this new parametrization A0A2A0 does not have 
the same effect on z p ,w p as A2A0A2 (nor do A1A2A1 and A2A1A2), so we no 
longer have the problem referred to at the beginning of this section. 

If we identify Riemann sheets that have the same values of z p ,w p for 
given t p , then Q ceases to be a Cayley tree and becomes the two-dimensional 
honeycomb lattice of Figure E3 

To see this, note that if Zp,Wp are the values of z p ,w p on the central 
sheet "D, then it follows from (|4.13j) that on any Riemann sheet the analytic 
continuations of z p , w p (for a given value of t p ) are 

z p = x m (z^ 1 , w p = x n (w ^ 1 , (4.14) 
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X [O] (2,2) 



1 (-2,0) 




2j (-3,-2) 



1 (-2,-3) 



Fi gure 3: The honeycomb lattice formed by the hypcrelliptic variables z, w. Circles (squares) 
denote sites of even (odd) parity. 



choosing the upper (lower) signs on sheets of even (odd) parity. Here m, n are 
integers satisfying 



m + n = (mod 3) on even sheets 
= 1 (mod 3) on odd sheets 



(4.15) 



The Riemann surface for z p , w p therefore corresponds to a two-dimensional 
graph Q, each site of Q being specified by the two integers m,n and corre- 
sponding to a Riemann sheet of the surface. 

This Q is indeed the honeycomb lattice shown in Figure 01 Adjacent sites 
correspond to adjacent Riemann sheets. Sheets of even parity correspond to 
sites represented by circles, those of odd parity are represented by squares. If 
i is the integer shown inside the circle or square, then on even sites x p ~ u> 1 , 
and on odd sites y p ~ u % . The numbers shown in brackets alongside each site 
are the integers m,n of l|4.14|) : we refer to the corresponding sheet as "the 
sheet (m,n)". 

We can trace this reduction of Q to the fact that the automorphisms Ai 
applied to z p ,w p ,pj satisfy 



AiAjAk = A^AjAi for all i,j,k . 



(4.16) 



For instance, the relation ^^l^o = AqAiA2 implies that D221 = Don, and 
indeed we see from Figure 01 that these are the two three-step routes from T> 
to Y. Similarly for the two routes to Z, or the two routes to X. 
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We noted in the previous section that the function T2(p,q), considered as 
a function of either p or q, also satisfies (j4.16|) . It follows that its Riemann 
surface can be embedded in that of the reduced graph Q. 

We focus on the q dependence of T2(p,q). Let its value when p,q £ T> 
(given by eq. 12. 7j) be r 2 {p, q)- Then by iterating IjH.fij) we find that on the 
Riemann sheet (m, n) 



T 2 (p,q) = v; g T°(p,q) , (4.17) 



where 



r = (m — 2n)/3 on even sheets, 

r = (m— 2n — 1)/3 on odd sheets. (4-18) 

(c) The function T(p, q) 

We can eliminate the distinction between even and odd sheets by working, 
not with T2 (p, q) , but instead with the closely related function 

TM = 44 = . T / P '"f v2 , , (4.19) 
a{p,q) T2{p,Vq)T 2 (p,V 2 q) 

(using eqn. 12.10(1 . Then (j.3.6j) becomes 

T(p,A t q) = v^W-^Tfaq) (4.20) 

for i = 0, 1, 2. If Tq(p, q) is the value of T(p, q) in V, then it follows that for 
a given value of t q , 

T(p,q) = v™- 2n T (p,q) , (4.21) 

on all sheets (m,n). 

It follows that T(p, q) is a single- valued meromorphic function on the Rie- 
mann surface G, and that the orders of its zeros and poles are linear in the 
integers m, n of (|4.14|l (with p replaced by q therein). This suggests that it 
may be possible to write T(p, q) as a product of functions of z q ,w q , and indeed 
we shall find that this is the case. 

We shall explicitly exhibit the dependence of T(p,q) on z q ,w q by writing 
it as T p (z q ,w q ). Then the three relations (|4.20|) become 

T p (x~ 1 /z q ,X~ l /w q ) = Tp(x/Z q ,l/W q ) = V pq T p (z q ,W q ) , 

T p (l/z q ,x/w q ) = Vp q 2 T p (z q ,w q ) . (4.22) 

If p is fixed within P, then from H2.7|) and (|4.19|) . T p (z q ,w q ) is bounded 
and has no zeros or poles for q € T>. (Both r| and a become infinite as 
Uqi~kq ~^ °°> but their ratio remains finite.) Together with this restriction, the 
relations (|4.22j) define T p (z q , w q ) to within a multiplicative factor independent 
of q. To see this, suppose we had two such solutions, then Q4.22JI would imply 
that their ratio was unchanged by q — > Aiq, for i = 0, 1, 2. This would mean 
that the ratio, considered as a functions of t q , did not have the branch cuts 
Bo, B±, £>2- It would therefore be an entire bounded function of t q , and hence 



11 



by Liouville's theorem a constant (independent of q). This constant could be 
determined from the product relation (12.10(1 , which now takes the simple form 

Tp(z q ,Wq)T p (-l/wq,z q /w q )T p (-Wq/z q ,-l/z q ) = 1 . (4.23) 

We shall also use the p-relation (|3.9|) . Together with (|3.1Uj) this implies 
that 

T(p,q)T(A iP , q ) = & - 6 . (4.24) 

5. Various p, q relations 

First we present various relations that enable one to express certain ratio- 
nal functions (including v pg ) of x p , . . . ,t q as products of elliptic functions of 

Z p , W p , Zq , W q . 

Some can be obtained by applying the rule (|4.ip to eqn 34 of (B 1993a) 
or to eqn 4.9 of (B 1993b). In particular, we obtain 

^ I gfej = <KV%)0KM>) • (5-1) 

Applying the automorphism U to the q variable in (|5.1j) . and then using 
the V and R automorphisms, we obtain the two sets of relations 

fi p (y q - UJ j X p ) _ j-i(t>{x m Pm-i+l Qm-j-l) 



Hq(y p - V l X q ) (f>(x m p m -i-iq m - j + l) 

HpHqjXg - Ujl Xp) = ^ m _j (/>(a3 + 1 p m _ i+ i/g m _j_i) 
yp-ujiyq 4>(x3- 1 p m ^ i - 1 /q m -j +1 ) 

true for all integers m. 

The function (f>(z) has a leading factor z 1 / 3 , which gives a contribution 
z q Wq/(z p Wp) 1 ^ 3 to the RHS of (|5.1|) . This cube root can and should be chosen 
to be 93/^3. The corresponding contributions in (|5.2() are p m -i/q m -j and 

Pm—i Qm—j ■ 

We define the elliptic function 

00 

g r (z) = Y[(l-x rn - r z)(l-x rn /z) (5.3) 

71=1 

for integer r (in particular r = 1 and 3), satisfying 

g r {x r z) = g r (l/z) = -z~ 1 g r (z) . 

Then h(z) = uj 2 g 1 (uj z) / 'gi{uj 2 z) and <j) h {z) = g 3 (x/z)/g 3 (xz), so (jP)) . (jnU|) 
can be used to relate various gi and 53 functions. In particular, from (|4.8|) 
and (UTTTl) . 

93(pj)/g3{p~i) = ^plg^p^/g^Pi) , (5.4) 

for any cyclic permutation k} of {1,2,3}. We note that (|4.2j) can be 
written 

- k' 2 = 27x/g 3 (x) 12 (5.5) 
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cinci set 

7 = (-27xA; 4 A ,2 ) 1/6 = -3uk 2 / 3 / gi (u) 2 

= k 2/3 g 3 (x) 2 = l + 7x + 8x 2 + 22x 3 + 42x 4 + --- . (5.6) 

It follows from ()5.3|) and (|5.5|) that these various expressions for 7 are all 
equal. 

We also define 



, v (y 3 p - x 3 q )(y 3 q - Vyg 

Vl[P ' q) (y p 3 - y 3 )(x 3_ 1/y 3) 



and note that 



= vi(p,q)v 3 (p,q) . (5. 



We find the following six identities and indicate below the method of their 
proof. The first two can also be derived from Q5.2[) by applying the dual- 
ity/conjugate modulus mapping of the Appendix. 

We have also checked all the identities of this section numerically, for 
arbitrarily chosen x, z p , z q , to 25 digits of accuracy. 

y q - tn j x p fipiJ, q _ -_j h(uj m p m -i-iq m - j+1 ) 



y p - ^XqflpHq h(uJ m p m - i+ iq m -j-i) 
X q H q ~ U j XpH p ^ /l(^ _1 p m -i_l/g m _j +1 ) 



y p fj, q - ujtyqfip h(ujj +1 Pm-i+i/qm-j-l, 



(5.9) 



M P ,^U^y (5-10) 

f X Z l W \ TT 93(-X 2 ZqWqPj) 

V3{P,q) = - s — 7 jztt- , 5.11 



u 2 gi{pi qj) gi(qj/Pi) 93(qj/Pi) 9z{x 2 Pi Qj) 



gi(^ 2 Pi)9i(^ 2 /Pi)9i(^ 2 Qj)gi(^qj) gz{~\-IPi)g?,(x/p i )g 3 {q :j )g 3 {x 2 q j ) 

(5.12) 



j (u;% - ujH q )(cu l+ H p t q - 1) 



for i, j, m = 1,2, 3. 

Letting q — > Vg, we see that (|5.1()|) remains true if z q therein is replaced 
by —1/wq or by —w q /z q . Also, (j5. is true if —z q w q is replaced by z q /w 2 or 

q/ z q- 



by -w q /z 2 



a) Proof of the identities 
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We have proved the identities (|5.9|) - (|5.12|) . Let E q be the ratio of the 
LHS to the RHS of any of these identities. 

First note that a priori we expect E q to be a single-valued function of t q 
only if we introduce the branch cuts Bi of Fig. ^hito the complex i g -plane. 

Indeed, the identities l|5,9j) involve \x q and qj, and these will need additional 
cuts linking the Bi in order to completely fix the choice of the cube root in 

(HE}. 

The map that takes q from one such choice to another is the map q — > Mq. 
However, from (|2.3j) and (|4.12|) . this merely increases (decreases) i,j,m by 
one in the first (second) set of the identities ()5.9|) . If we take some symmetric 
function & q of all the E q (for all i,j,m), then $ g will be invariant under 
q — > Mq, which means that these extra cuts are not needed for this function. 

The next step is to use ()3.2[) and (|4. 13f) to show that the mapping q — > Aiq 
merely permutes, and possibly inverts, the E q of each set of identities. This 
is true for each of the sets of identities (|5.1UI ) - ()5.12|) , but for (|5.9|) the two 
sets are interchanged by this mapping. 

Again, let & q be some symmetric function of all the E q (and if necessary 
their inverses) within a set (e.g. the sum of the fifth powers of every E q in one 
of the eqns. 15.121 summed over i and j), now regarding the two identities (|5.9|) 
as forming one combined set. Then $ g will be unchanged by each of the three 
mappings q — > Aiq, for i = 1,2,3. This means that it has the same value on 
either side of any of the branch cuts Bi. The cuts are therefore unnecessary: 
<5 g is a single-valued function of the complex variable t q . 

The only possible singularities of <fr g are therefore poles, arising from poles, 
and possibly zeros, in the E q . The only places these can occur are when 
tg = tp, t q = 1/tp, t q = and t q = oo. We can restrict our attention to 
p,q € T>, since Q q is unchanged (for a given value of t q ) by crossing the Bi. 
Thus \n q \ and \n p \ are both greater than one, and x q ,x p each lie in a region 
near the points x q = 1, x p = 1. This means that the only points to consider 
are: 

1) {x q ,y q ,Vq,qj,qj,qj} 

2) {x q ,y q ,H q ,qj,qj,qj} 

3) {xq,y q ,Vq,qj,qj,qj} 
^){xq,y q ,Vq,qj,qj,qj} 

for all i and m. 

Each E q can be written as a ratio of pole-free functions. In every case, if 
the numerator (denominator) has a zero at one of the points (|5.13|) . then so 
does the denominator (numerator), and both zeros are simple. Thus no E q 
has a pole or zero at any of the above points. 

The function $ g is therefore a single-valued and analytic function in the 
complete t q plane, including the point at infinity (this is the last of the points 
listed above). By Liouville's theorem it is therefore a constant. 

One can write down a polynomial of finite degree whose roots are the E q 
and whose coefficients are symmetric functions Q q . Since each such coefficient 



= {x q ,uj-y q ,u) fi q ,q j+i ,q j+i ,uj q j+i \ , 
= {l/x q ,u l /y q , uj m ~ 1/2 x q fiq/y q , 

qlU'V-i-i'^^-i-i} > ( 5 - 13 ) 
= {k 1 / 3 ,0,oo m /k >1/3 ,u J 2 ,l,Lo m x 1 / 3 } , 
= {A; 1 / 3 , 00, w m (— fc/fc') 1 / 3 ^, w, 1, uo m ~^x 1 / 3 } , 
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is a constant, so are the roots. Thus every E q is a constant. By looking at 
special values of q, e.g. one of the values above, in each case we can show that 
the constant is unity, and hence prove the identities - l|5.12|) . In the last 
two identities one needs to take the limit as y q — > or oo, using (|4.2() . (|4.9|) . 
jOffll and (|52J - CT . 

(One can streamline the procedure: for instance, each of the 27 identities 
in the first set of relations (|5.9ft can be obtained from one of them by using 
the mappings q — > Vq, q — ► Mg, p — ► Vp, p — > Mp, so it sufficient to do the 
last step for just one of the equations, say i = j = m = 0.) 



6. Calculation of T(p, g) 

We now look for solutions of (j4.22j) , using the identities - (|OT|l . This 
leads us to define, similarly to section 6 of (B 2006), the function 

°° (I - r r i z V 

FJz) = TT ,) . (6.1) 

It satisfies 

F r (z) = l/F^z" 1 ) = -z- 1 g r (z)F r (x r z) (6.2) 

and is a natural extension of the elliptic function g r {z). 
We further define the functions 

3 

R x {z) = l[[F 1 (z/p j )/F 1 (zp j ) , 

3 

suppressing their dependence on p. They have been constructed so that 
Rl(x~ m /q l )=v 1 (p,q) m R l (q i ) , 

Rz{x l - Zm /qd = (-fc) m x^ 3 ™" 1 )/ 2 v 3 (p, qrSafa) , (6.3) 

for i = 1,2,3 and all integers m. Because they are products over pi,P2,P3 
(or p 1} p 2 , p 3 ), they, like vi(p, q), v 3 (p, q), T(p, q), are unchanged by p -> Fp. 
They only have zeros or poles when vi(p,q),V3(p,q) have zeros or poles, i.e. 
when x 3 or y 3 equals x^ 3 or y^ 3 . None of these zeros or poles occur when 

P ,qev. 

The factors (—qi) m , x m( - 3m ~ 1 ^ 2 in (|6.3|) are irritating as they do not occur 
in Q4.22JI . However, they are independent of p, so we may hope to remove them 
by introducing some additional simple factor Xq that is also independent of p. 
Also, (|4.22j) is unchanged by multiplying T(p, q) by any function r\ v of p only. 

We therefore try the ansatz 

3 

T p (z q ,w q ) = ^^n^fer^fe) 6 ' , (6.4) 
i=i 
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where a\, a 2 , a 3 , 61, b 2 , 63 are integers, to be determined. We first seek to 
satisfy the relations (|4.22|) . and find that we can match the powers of v pq = 
vi(p,q)v 3 (p,q) therein by taking 

0,2 = ax — 1 , a 3 = ai + 1 , b 2 = b 3 = 1 , 61 = -2 . (6.5) 

At this stage we are free to choose 01, which corresponds to multiplying 
T p (z q ,Wg) by a factor R 1 (q 1 )R 1 (q 2 )R 1 (q 3 ) = R\{z q )R 1 (-l/wq)Ri{-w q / z q ). 
This factor has no effect on the relations (|4.22() and is bounded, with no zeros 
or poles, for p,q G D. From the argument after (|4.22() . it must therefore be 
independent of q. Taking t q = 00 and qi = to, or t q = and qi = cu 2 , we obtain 
the identity 

3 

YlRliqj) = Ri(u>) 3 = Ri(uj 2 ) 3 (6.6) 
i=i 

for all q. 

Without loss of generality we can therefore choose a\ = 0, giving a 2 = 
— 1, a 3 = 1. Substituting the ansatz (|6,4j) into (|4.23j) . the R\ and R3 functions 
cancel out, leaving 

Vp Xq XVq XV 2 q = 1 ■ ( 6 - 7 ) 

Since p and q are independent variables, r/ p must be a constant (independent 
of both p and q). We can absorb this constant into the factor Xq i n <|6-4J) . so 
we can set 

r, P = l . (6.8) 

We can calculate Xq from (|4.24|) . Take i = 1 therein, so = A\p and 
= x 2 ~ j /pj and = x 1-315 ^' 2 ) /pj for j = 1, 2, 3. Let be the function 

R m (z) defined above, but with p replaced by p'. From the above definitions 
and properties we can verify that 

- pm9l{zpMz/pi) 



gi{zp3)gi(z/p 3 ) 
and 

R' 3 (v)R 3 (v) = 1 _ . (6.9) 

gz[v/p 2 )g 3 (x i vp 2 ) 

Using these, we find (after some work) from (|5.12j) and (|6.4[) that (|4.24j) is 
satisfied for i = 1 iff 



2 g3{q 2 )g3(x q 2 )g2,{q 3 )g 3 {x g 3 ) 
Xq gziqiYg^x 2 ^) 2 



Taking i = 2 or i = 3 in (|4.24() merely permutes Pi,P2,P3 in the working, 
which leaves the p-independent result (|6. 1U|) unchanged. 

We can use (|4.1Uj) to eliminate the factors gs(qj) in favour of gs(x 2 qj). 
Using also ()6.7j) . we obtain 

_ g M^%) 

q 1 g 3 {x 2 q i y 



1(3 



Writing Xq a s x{ z q-, w q) an d using (|4.1U|) . we find that 



-xq 3 x(x 1 /z q ,x 1 /w q ) = -xq 2 x(x/z q ,l/w q ) 

= X~ 2 qi 2 X{t/Zq,x/w q ) = x{Zq,W q ) . (6.12) 

Our ansatz (|6,4|) is now 

rr, ( \ -Rl(g3)^3(g 2 )^3(g 3 ) , R io\ 

Tp{Zq ' Wq) = Xq RMRM* ■ (fU3) 

Using (|6.3|) and (|6.f 2|) . we find that this expression does indeed satisfy the 
relations (|4.22|) , so from the argument following (|4.22j) we see that the expres- 
sion 1)6.13(1 must be correct to within a factor that is independent of q. It also 
satisfies (|4,23|) . so this factor must be a cube root of unity. Since (|4.24|) is 
satisfied, this root must be unity itself. This completes the proof of 1|6.13|) . 

(a) Relation to the order parameter function S(t q ) 

In an earlier paper (B 2006), the author considered the function S(t q ) that 
occurs in the derivation of the order parameter of the chiral Potts model. This 
is given by 

log S(t q ) = -^ogk + — j o log [A(0) - t q ] d9 , (6.14) 

A(8) being defined by (|2.8|) . This is very similar to the function T2(p,q) of 
this paper, in fact if qz is the point (3) of eqn. 1)5. 13|) . with ^ = 1/k' and 
t q = 0, and Q4 is the point (4), with n q = t q = oo, then 

s[tq) = kW-nM ■ (6 - 15) 

For N = 3, we can use our results (|4.19|) . (|6.13|) to express the RHS of 
(|6.15|) in terms of our elliptic-type functions. Because qj = 1 for both points 
(3) and (4), the R3 factors cancel. The R\ factors reinforce, giving 



where 

G{z) = F 1 (ujz)/F 1 (uj 2 z) . 

This is indeed the result given in equations (57), (62) of (B 2006) (with p 
therein replaced by q). Because the R3 functions have cancelled, it contains 
no elliptic-type functions with arguments z q /Wq, —w q /z q or —z q w q . 
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7. Summary 



For N = 3, we have written the integral expression (|2.7j) as a product of 
generalized elliptic functions, with arguments that are ratios and products of 
the variables pj , qj , pj , defined in l)4.5|) , Q4.7JI . These are the variables of the 
alternative hyperelliptic parametrization of the chiral Potts model. 

Is this progress? To the author it seems that the answer is indeed yes: 
functions such as the F r (z) of 1)6. 1JI occur naturally in the free energy of other 
models, such as the Ising, six and eight-vertex models, so it is interesting 
to see them occurring again in the three-state chiral Potts model. Indeed, 
it may indicate some intriguing relations between such models. Certainly it 
provides an explicit formulation of the meromorphic structure of T2 (p, q) on 
its Riemann surface. 

Unfortunately it is not clear how one could proceed to N > 3. There seems 
then to be no reason to expect to be able to express quantities such as v pq as 
products of single-argument Jacobi elliptic functions. The best one can hope 
for is to write them as ratios of hyperelliptic theta functions, as in (B 1991b). 
These are entire functions of two or more related variables, all of whose zeros 
are simple. To write T%(p,q) in such a parametrization, one would need to 
generalize the hyperelliptic theta functions to make the order of the zeros 
increase linearly with the distance of the zero from some origin. The author 
knows of no such generalization. Indeed, while Jacobi's triple product identity 
enables one to write Jacobi theta functions as either products or sums, there 
seems no reason to expect the same of functions such as the numerators or 
denominators in 1[6.1|) . 

There is also a problem with extending our working to the free energy of 
the full N = 3 chiral Potts, which is given by the double integral in eq. 46 
of (B 1991a) and eq. 61 of (B 2003a). The author showed in (B 2003b) that 
the graph of the Riemann surface of this function has one more dimension 
than T2(p, q). If we fix p and consider the free energy as a function of q, then 
we do not have the relation (|4.16|) and we need a three-dimensional lattice to 
represent the Riemann surface, rather than the honeycomb lattice of Figure 
IHI This implies the need for one more "hyperelliptic" variable, in addition to 
z q and w q , and it quite unclear what this may be. (No such extra variable is 
needed for the N = 2 Ising case.) 

Appendix A: The duality map 

The mappings R, U, V, M leave k,k',x unchanged. There is another map 
that takes k' to 1/k' and p to p', where 

{ay>%>'>/V> V) = {w 1/4 aw w 1/4 y p //i p , 1/ /j, p ,uj 1/2 t p } . (Al) 

Let x = e~ nX , z = e ma and write the RHS of (jJ2J) as r(A). Also, write 
h(z),4>(z) as h(a, A), <f>(a, A). Set 

A' = 4vr/(3A) . (A2) 



18 



Then, as in §15.7 of (B 1980), one can establish the conjugate modulus 
relations 

h(a,X) = 0(-2ia/A,A') , 4>(a, A) = /i(2ia/3A, A') . (A3) 

Also, r(A) = l/r(A'), so replacing A by A' does indeed invert k', as given by 
(|4.2|) . Further, if pj = exp(i7raj), then we can choose 

a.\ + ct2 + «3 = — 2 . 

(This enables one to take ctj = —2/3, for all j, for the interesting case when 
y v = and pj = to^ 1 .) 
If we set 

2i 

«j = -ttj-i , (A4) 

then it is also true that a[ + a' 2 + a' 3 = — 2. In fact the mapping A, ctj — > A', ct^ 
is self-reciprocal. Further, using (jA3|) . we find 

h(otj, A) = 4>(a'j + i — a'j_ 1 + i\' , A') , (j)(aj+i — a,j-\-\- iA, A) = h(a'j, A ) . 

The left-hand sides of these last two equations are the functions h(pj), 
4>(xpj) of eqns. (|4.9jl . (|4.10|) . It follows that these two equations are inter- 
changed by the duality mapping (|A1|) . provided we also use the conjugate 
modulus form (jA3|) of the elliptic functions and then make the transforma- 
tion A, aj <-> A', ctj. 
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